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SUMMARY 

Analysis of rotating stall of coii 5 )ressar Blade rows requires speci- 
fication of a dynamic lift curve for the airfoil section at or near 
stall, presumahly including the effect of lift hysteresis. Considera- 
tion of the Magnus lift of a rotating cylinder suggests performing an 
unsteady boundary- layer calctilation to find the movement of the separa- 
tion points of an airfoil fixed in a stream of variable incidence. Then 
consideration of the shedding of vorticity into the wake should yield 
an estimate of lift increment proportional to time rate of change of 
angle of attack. This increment is the amplitude of the hysteresis 
loop. 


An approximate analysis is carried out according to the foregoing 
ideas for a 6:1 elliptic airfoil at the angle of attack for maximum 
lift. The asB UTnp ti nn of small perturbations from maximum lift is made, 
permitting neglect of distributed vorticity in the wake. The calculated 
hysteresis loop is counterclockwise. The computed increment of lift 
coefficient is quite large, indicating appreciable unsteady lift hyster- 
esis for quite small reduced frequency of flow oscillation. It is 
assumed that to the order of this analysis, the wake begins at the 
separation point defined by zero shear. This assuD 5 )tion is questionable 
for unsteady flow. 

Finally, a discussion of the forms of hysteresis loops is presented; 
and, for small reduced frequency of oscillation, it is concluded that 
the concept of a viscous "tirK lag" is appropriate only for harmonic 
variation of angle of attack with time at mean conditions other than 
maximum lift . 



2 


MCA TN 3571 


nmiODucTiON 

Tlie phenomena of "stall flutter" and "rotating stall" which may 
appear in an axial-flow cor^ressor both involve fluctuations in flow 
about blades operating near their aerodynamic stall point; that is, at 
an average flow Incidence angle near that corresponding to maximum blade 
lift. 


The analysis of stall flutter has been held back by uncertainty as 
to the dependence of airfoil lift and moment on a fluctuating incidence 
angle near stall. Msndelson in reference 1 proposed assun^jtion of the 
linear aerodynamic force and moment relations appropriate to steady 
flow at small angle of incidence, 2 nodlfied by the further assumption 
that, as the airfoil oscillates, the forces and moments lag behind the 
angular displacement of the airfoil, owing to viscous effects. Such a 
time lag represents an unsteady hysteresis which may provide cyclic 
work to amplify or maintain flutter. Perhaps the first experimental 
study of lift hysteresis was that of Farren (ref. 2). Half man, Johnson, 
and Haley (ref. 3) and Schnlttger (ref. 4) have more recently studied 
aerodynamic hysteresis experimentally and have presented empirical anal- 
yses of their results. 

An analysis of rotating stall (e.g., that of Sears, ref. 5) also 
requires specification of a dynamic lift-incidence relation (or the 
equivalent, as in the study of Emmons, Pearson, and Grant, ref. 6, and 
in Marble's analysis, ref. 7). Sears has adopted Mendelson's phase-lag 
hypothesis, and this phase angle is an undetermined parameter of his 
analysis . 

The concept of viscous time lag is not entirely satisfactory, 
however, partly because the phenomenon itself is unexplained, but chiefly 
because the concept obviously cannot describe a lift-hysteresis loop 
which might occur at a nominal condition of maximum lift. 

The phenomenon of aerodynamic hysteresis presumably depends, at 
least in part, on the airfoil boundary layer. Also, in this study, 
hysteresis is taken to be a fundamentally unsteady phenomsnon, not ex- 
plainable by consideration of the steady or quasi-steady boundary layer. ^ 
In the present report, consideration is given to the incompressible 
flow field about a single airfoil fixed in a flow of oscillating 


^Airfoils with lift curves which break sharply at stall may show 
lift hysteresis in steady flow, a phenomenon distinct from that under 
study herein. 
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incidence^ under thB aBsun^jtion of an unsteady but nearly quasi-steady^ 
laminar boundary layer. (This sort of boundary layer is analyzed in 
ref. 8.) The analysis of this flow field is undertaken in order to gain 
an understanding of the cause of lift hysteresis and to describe its 
form and (crudely) its magnitude for a special airfoil at maximum lift. 

The basic quasi-steady flow to be used in the present analysis is 
provided by Howarth's analysis (ref. 9) of the way the laminar boundary 
layer (and hence circulation) about an infinite elliptic cylinder de- 
^ pends on angle of attaok, applied at a condition of maximum lift. In 

reference 9, Howarth makes assun^tions under which his result becomes 
quantitatively inaccurate^ though both the approach and the result are 
qualitatively instructive. The same limitations affect the present 
analysis . 

Holding the airfoil fixed while the flow direction oscillates sim- 
ulates the accepted picture of rotating stall, in which successive blade 
passages stall progressively along a perfectly rigid cascade. The some- 
what different case of an oscillating airfoil in a uniform stream, which 
is appropriate to the stall flutter problem, is not analyzed in this 
report. However, there is an example in the "oscillating airfoil" cate- 
gory which Illustrates the considerations underlying the present study; 
namely, the rotating circular cylinder in a uniform stream. If a circu- 
lar cylinder is fixed in a uniform stream, it of course experiences no 
lift. Further, if it is given an angular displacement, its lift does not 
change, but remains zero. Thus, this degenerate "airfoil" may be said to 
be in a stall condition, at maximum lift, in fact. Now, if the circular 
cylinder is given a constant angular velocity of rotation about its axis, 
then a circulation develops and an aerodynamic force (Magnus force) 
transverse to the flow direction is exerted. If the stream velocity is 
from left to right and the rotation is clockyise, then the force is up- 
ward (lift) . If the rotation is counterclockwise, the force is 
downward. 

This phenonenon is explained (ref. 10, par. 27) by consideration of 
the boundary layer. In the case of clockwise rotation, the upper surface 
of the cylinder is moving with, and the bottom surface against, the flow. 
Consequently, if circulation remains zero, the boundary layer separates 
later on the top and sooner on the bottom than is the case when the cyl- 
inder is not rotating. On the top, later separation means that the 
velocity outside the boundary layer is lower at separation. Now, the 
separation point signifies the beginning of a wake. Therefore, the 


O 

In a quasi-steady flow, quantities vary slowly' enough so that 
steady-state results apply at each instant of time although slight 
variations are permitted from one instant to the next. 
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clockwise vorticity shed into the wake, being proportional to the local 
outer velocity, is less on the top, and the counterclockwise vorticity 
shed at the bottom separation is greater, than in the case of no 
rotation. 

Therefore, owing to clockwise rotation, a net increase of counter- 
clockwise vorticity is shed- By the law of conservation of circulation, 
the circulation therefore cannot be zero, and a clockwise circulation 
must develop about the airfoil to con 5 )ensate for the shed vorticity. 
According to classical hydrodynamics, this circulation results in lift. 

If, instead of rotating steadily, the cylinder undergoes a rota- 
tional oscillation, the same considerations apply, if the reduced fre- 
quency of oscillation is small. In that case, the oscillating lift is 
proportional to the instantaneous velocity of rotation. Thus, when the 
"angle of attack" of the cylinder is increasing, there is positive lift, 
and when the angle is decreasing, there is negative lift; over a complete 
cycle, the curve of lift against angle of attack would be a loop. 

Therefore, the circular cylinder undergoing a rotational oscilla- 
tion exhibits lift hysteresis, by reason of the response of the boundary 
layer to the movement of the surface. In the more complicated problem 
of a noncircular cylinder, or airfoil, similar considerations may be 
expected to apply. Of coTxrse, in the airfoil problem conten 5 )lated in 
the present study, the acceleration of the flow field may be expected to 
provide an additional coiiponent of pressure lift, derivable from con- 
sideration of Kelvin's impxilse. 


PEEUMINARY CONSrDERATIQNS 
Statement of Problem 

Consideration is given to the lift of an isolated airfoil in the 
form of an infinite elliptic cylinder with a semichord 1 and a thick- 
ness ratio P, at a stalling angle of attack a to a stream of velocity 
U (see fig. l) . A full list of notation is provided in the appendix. 

Holding the airfoil position and the magnitude of U fixed, the 
angle of attack a is permitted to vary with time. Such a flow may be 
constructed by allowing a moving source Q to approach an airfoil fixed 


3778 



3778 


KA.CA M 3571 


5 


in an ottierwiBe uniform stream, the direction of approach being normal 
to the stream direction, as in sketch (a) : 



■jV'i i- 


Sketch (a) - Airfoil and source 


So long as h, the instantaneous distance from Q to the airfoil, is 
much greater than I, the airfoil finds itself effectively in a uniform 
stream of incidence 



+ 


Q 

2rthU 


and. of magnitude differing from U only to second order in Q/2jthU. 
The rate of change of angle of attack is 


a 


Q 


(-h) 


where the dot signifies differentiation with respect to time. The fore- 
going model applies qualitatively to the phenomenon of rotating compres- 
sor staJJ., if the moving source Q is taken to represent the approach 
of a flow blockage propagating along a cascade. 

The present analysis will be ceirried out as though a is a small 
constant. Actually, if a is quite small, and higher derivatives such 
as a are negligibly s mall , the analysis will be correct at each instant 
using the appropriate instantaneous value of a. (This is the first 
refinement over the quasi-steady assumption which uses instantaneous 
values of a itself; or for oscillatory a, the linear term of a Taylor 
series in reduced frequency.) It is clear that the solution of the prob- 
lem to order a provides a measure of hysteresis: Suppose an expression 

for lift is obtained in the form 

Cl = Cj + AuCj + oCj . + • • * 
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The first tvo terms are the qtiasi-steahy contrihutions . The term in a 
provides that, if a is in the process of increasing with time, then 
the lift is higher (assuming positive) than the quasi-steady value. 

The converse is true if a is decreasing. Thus, if a performs an 
harmonic oscillation, the lift curve is a loop lying to either sid^ of 
the quasi-steady curve of width proportional to the frequency of 
oscillation. 

In order that the present analysis hear on the question of rotating 
stall, nominal angle of attank must he selected for which the airfoil 
is in a stalled condition. Maximum lift is the most sin 5 )ly described 
stall condition. Accordingly, the nominal angle of attack is chosen as 
that for which the lift is a maxiimim. This selection is made for two 
more compelling reasons: 

(1) The result will tend to isolate the effect of hysteresis, in- 
asmuch as no quasi-steady change in lift results from change in a about 
the maximum lift value. Of course, if lift hysteresis is found under 

a mean condition of maximum 11 f t, then the idea of viscous time lag will 
thereby be shown to be inappropriate. 

(2) Any other assumption would lead to great theoretical complica- 
tion. The analysis is to be a perturbation of quasi-steady flow. If, 
at the nominal angle of attack, change in a resulted in a quasi-steady 
change in circulation, then, to the order of the present analysis, 
induced wake effects would require consideration. 


Potential Flow 

Outside the boundary layer of the ellipse, irrotatlonal incompress- 
ible flow is assumed. At each instant,^ the velocity potential on the 
surface of the ellipse is (ref. 11, par. 7l) 

(p = 7U(l + 3 ) cos (q-a,) - (l) 

where the surface is defined by 

X = 2 cos qj y = 32 sin q ( 2 ) 

Along the surface, measuring s clockwise, 

_ Bcp dq 
1 ” 5s^ ~ 5^ ds 

^ = - (sin^q + 3 ^cos 2 q)"^/^= _ 1 ( 3 ) 


^The quasi-steady assumiptlon applies precisely for the calcixlatlon 
of the unsteady velocity potential. 
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3778 


KACA TN 3571 


7 


At this juncture, the condition of maximum lift has not "been Imposed, 
and the circulation r is left unspecified. Of course, "both a and 
r may vary with time. 

The foregoing description of the potential velocity distribution 
is made on the assun^tlon that the boundary layer is negligibly thin 
everywhere on the ellipse. This assumption is usually quite proper 
ahead of the separation points. The assurption that the wake aft of 
separation does not luportantly affect the potential flow is not proper; 
certainly, this assumption is quantitatively poor at maximum lift, es- 
pecially if the potential flow is used to compute separation point. 
However, the results obtained on the basis of this assumption are ex- 
pected to have qualitative validity. 


BOUKDAEY-LAYER AHALYSIS 
Quasi-Steady Boundary Layer 

As a basis for subsequent calculation of unsteady effects, the 
quasi-steady laminar boundary layer on the ellipse may be approximately 
determined by the Karman-Pohlharisen integral method, as improved by 
Holstein and Bohlen (see Schllchting, ref. 12, ch. XIl). The differen- 
tial equation is 


dZ F(x) 
ds “ q ^ 



(5a) 


subject to the initial condition at the stagnation point (q = O) : 


Xq = 0.0770 (5b) 

where Z = 9 being the momentum thickness. The function F(x) 

is tabulated in reference 12. 


Determination of Z and r e.t maximum lift . - Given the velocity 
distribution q, the growth of the boundary layer may be computed from 
equations (5), the calculation proceeding until both the separation 
points (¥ and £, fig. l) are reached, for which (x)g = (x)-g- = -0.1567, 

In the present problem the potential velocity distribution has not 
yet been fully prescribed, since T remains unknown. FoUoid.ng Howarth 
(ref. 9), suppose that, for a given a, a value of r is assumed, and 
the boundary-layer calculation is carried out. Then, at the top separa- 
tion point, clockwise vorticity is shed into the wake at the rate 


u 


Su 


1 r 2 , 
= 2 


i6 


( 6 ) 
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wtiile at 8, counterclockwise worticity is shed at the rate ) . 

Now, if (u]^)g differs from (u]_)-g-, there is a net change of circulation 

in the wake, which is not possible in steady flow. Therefore, new trial 
values of F must he assumed, until the particular value of F is de- 
termined for which the solution of equations (s) yields separation points 
for which (u^)g = (ui)-* 

The foregoing procedure must he carried out for each of a number 
of angles of attack in order to determine the values of a and F at 
maximum lift. In reference 9, Howarth performed these calculations for 
3 = l/S and determined a theoretical lift curve. The maximum value of 
F/2jttn was found to he 0.0761 at an angle of attack a = 7°. His com- 
plete dlstrihution of Z at this condition is not presented in ref- 
erence 9. Accordingly, the calculation has heen repeated, hy integrating 
equation (5a) in the form 

^=A(ti)f(x); x=a(Ti)Z (7a) 

where 

(7h) 


( 8 ) 

and the In- 


(7c) 

The solution was obtained using a step-hy-step method in which a 
parabola is passed through two known values and the next unknown value 
of dZ/dT), Integrating to find Z in terms of the unknown dz/di], then 
applying equation (7a) at the unknown point to solve for dZ/di]. Two 
starting values wege found from a Taylor series about tiq. The step 
size in t] was 10 except near the stagnation and separation point where 
finer spacing was used. The solution of Z is shown in figure 2 and in 
table I. The separation points, for which x = -0.1567, occurrfed at 
Ti = 80.0° and 340.83°. Of course, q should be the same at s’ and £. 
The difference cited in the table indicates the degree of error present 
in the calculations. 

Determination of 5z/^a . - For subsequent use in the unsteady equa- 
tions, it is necessary to know the rate of change of Z with a in 
quasi-steady flow. At maximum lift, when a is changed, the quasi- 
steady boxindary layer changes, and the locations of the separation points 


A(n) = - f ^ 

and, from equation (4=), 

1 

At the forward stagnation point where q = 0, = 190.74°, 

itial condition is, from equations (5), (6), and (7a), 

Zq = 0.0770/aQ = 0.00407 


n 

“-Bin (ti-7°) + 0.0761 
6 


61 
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are changed. Of course, the velocities at £ and s must remain equal, 
because Sr/Sa = 0 at maximum lift, by definition. 

Differentiating equation (7a) yields 


5 

dll 


^ = B(ti)Z^ + b(Ti) 


(9a) 


where 


B = J1 F'(x) . 
q 

b= -A — F+ 


Zq 


(9b) 


aq 


R 




The initial condition for Z^ at tJq is determined by specifying that 
dZ /dT) must be finite there; from equations (9a) and (9b), 






V' 


= -0.00660 


'0 


(9c) 


In equations (9) all qxaantities are to be evaluated at the condition of 
maximum lift; the appropriate superscript (o) is omitted for brevity. 

Of coxirse, for purposes of finding q^ and q^^^ in equations ( 9 ), the 
angle of attack of 7° should be replaced by a in equation (8) and set 
equal to 7° again, subsequent to differentiation. Equations (9) have 
been integrated to yield Zq_, by the same method as described for find- 
ing Z, and the result is shown in figure 3 and in table I. 


Unsteady Boundary Layer 

The next step in the analysis is to, determine the dependence of the 
boundary layer on the angular velocity a, assumed small. To this order 
of approximation, equation (s) describing the potential flow must be 
modified to Include the possibility of a contribution to circulation in 
proportion to a (or, in dimensionless form, 6 = aj/u), as follows: 

I sin(ii-7°) + 0.0761 +. (lO) 

The coefficient of circulation hysteresis j nmst be found from a con- 
dition of vortlcity shedding at the separation points of the unsteady 
boundary layer. The contribution to lift proportional to a then fol- 
lows. Determination of the proper vorticity-sheddlng condition will be 
discussed in a subsequent section. 
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The unsteady form of the K^man momentum equation (see ref. 12) is 




( 11 ) 


steady equations (5a) and (7a) are obtained from equation (ll) hy omit- 
ting the term on the right side. An iteration procedure might then he 
adopted: The quasi-steady S'*" can he substituted into the right side 

of equation (ll), and a new solution obtained, to include first-order 
unsteady effects. Using the definitions of reference 12, and writing 
the time derivative in equation (ll) as 


5t 


(u^6*) = a, (u^&*) 


( 12 ) 


there is found, corresponding to equation (5a) , 




ds 


q i‘ 


2 Z 


(13) 


The function f^('^) is tabulated in reference 12, and q is given by 
equation (lO) . Again, for purposes of finding q^^ and q^ , the angle 
7° in equation (lO) should ten^jorarily be replac^ by a. 


Instead of the Indicated iteration, in the present study the equiv- 
alent procedTire is adopted of finding the coefficients of the expansion 

Z = z(®^ + doZ„ + 6Z- + • • • 

CL & 


The coefficients z(o) and have already been found (eqs. (7) and 
(9) and figs. 2 and 3). The derivative Zg remains to be found. Dif- 
ferentiating equation (13) and noting that = qgZ^^ + q^^Z, 


dZ 

e 

ds 


F'(x) 


% 

q 


Zg + F' 







In equation (14) and hereinafter, evaluation of quantities in the steady 
state at maximum lift is to be understood, and the superscript (o) is 
omitted for brevity. 

It is iB5)ossible to ensure a finite value of dZg/ds at the stag- 
nation point because of the second-order pole -2f^Zq^q^. The physical 
reason for this result is the fact that, at the stagnation point of a 
certain instant, the boundary-layer velocity profile will not vanish, 
as in steady flow, but rather will respond to the instantaneous accelera- 
tion 6 more pronptly than the outer potential flow, A profile of 
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magnitude 6 may thus he expected, vanishing in the outer stream as 
well as the wall as shown in sketch (h) : 


/ 



Sketch (h). Nose of airfoil 
Therefore, the definition of momentum thickness 


0 



shows that, if the velocity u in the boundary layer has a part pro- 
portional to Un and a part proportional to e , then the part of 0 
(and hence of Z = which is proportional to e must have a 

simple .pole in u^_. Actually, of course, the quantity 0 is inappro- 
priate for defining a thickness of a profile of the type shown in sketch 
(h), and the appearance of a pole in 0 simply indicates this lack of 
physical significance. 


The foregoing considerations suggest that a new variable W be 


defined to replace Z 


' 6 - 


W = qZ 


(15) 


Substituting equation (l5) into equation (l4) yields 
ds a 


Cl + P')w + Zf'qeo - f ^ ^ 

(16) 


O 
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In order that dW/ds he finite at the stagnation point, the two poles 
on the right side of equation (l6) mast cancel, yielding the initial 
condition 


where the numerical, value is obtained using equations ( 3 ), (7c), and (lO), 
and the tables of reference 2. Inasmuch as Qg = t/R (eq. (lO)), the 
function W may be split into two parts, as follows; 

W = X + rY (18) 

so that (changing to as independent variable) equations (I 6 ) and (l7) 
provide 


a + F')q, 


= -0.001015 


iX7) 


where 


and 


^ = C(ti)X + c(ti) 


C(n) = (1 + F’) ^ 


c(ti) = 


RZ 


q 2Z 



(l9a) 


(19b) 



Xq = -0.001015 

(19c) 


C(n)Y + d(Ti) 

(20a) 

where 

d(Ti) = 1 _ ZF> ^ 

(20b) 

and 

0 

II 

0 

(20c) 

Equations 

(l9b) and (20b) are evaluated using the solutions for 

Z and 


Z , the tables of reference 12, and equations ( 3 ) and (lO) . 

Solutions of equations (l9) and (20), obtained by the method used 
to find Z, are presented in figures 4 and 5 and table I. 
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DETERMIMTION OF LIFT 


Unsteady Balance of Shed Vorticity 

In order to determine the unsteady pressure lift, the coefficient 
Y (eq. (lO)) must he determined. In steady flow (see the discussion 
acconpanying eq. (6)), the circulation F was obtained hy req^liring 
that vorticity he shed in equal and opposite amounts at the two separa- 
tion points. 

In the present imsteady problem, the net rate of vorticity appear- 
ance in the wake must vanish, not only in the quasi-steady approximation, 
hut to order a as well, in view of the assumption that the quasi- 
steady circulation is maximum. By the classical theorem concerning con- 
stancy of circulation, any net rate of discharge of vorticity into the 
wake must he balanced by a rate of increase of circulation about the 
body. If the airfoil is nominally at maximum lift, then the circulation 
terms (eq. (lO)), to order e = aZ/U, are 

0.0761 + rt 

The rate of increase of this expression is at most of order a, and 
therefore there cannot be any net discharge of vorticity to order a. 

If the airfoil were not at maximum lift, then the expression for 
circulation would contain a term proportional to Aa, which would change 
at the rate <x, and would have to be balanced by a net rate of vorticity 
discharge of order a. In turn, this distribution of circulation in the 
wake would induce further modifications of the potential flow. Therefore, 
the assunption of maximum lift permits the neglect of the induced effects 
of distributed circulation in the wake. 

Movement of separation points . - In order to effect a balance (to 
order a) of vorticity shed at the separation points, the movements of 
the separation points must be taken into account. The position of the 
top separation point ¥ may be written 


n = 


(q)^ + Ax 



( 21 ) 


The coefficient 


is 


At separation (x = -0.1567),’ 


o 

obtained from the quasi-steady solution: 




14 


MCA TN 3571 


From equations (7), 

Bx 

5^ 


= -i(= 






') 


vhence, holding x fixed at -0.1567, 





(22a) 

(22b) 


(23) 


where Z- and Z may be obtained from table I. 
a 


The coefficient (Sri/^e)— comes from the unsteady solution: 


(il = f4tl 


From equations (7), 

■>^6 = " R (^Igq^ + Zq^g) 

which, upon substitution of equations (lO), (l5), and (l8), becomes 
Xc = - r (X + rY) ^ ^ sin 


R 


q 4 

^ 72R 


2^ 


(24) 


Equations (22a) and (24) thus provide that 


ey_ 

s 


^ % f % 35Z . o 

^ q ^ ^ =!■> 21 ] 


%1 *W*^ 


( 25 ) 


At thfe bottom stagnation point, equations (2l), (23), and (25) 
apply, with subscript s repleicing subscript s’. 
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Equation of Torticity shed at separation points . - Taking- into ac- 
count the motion of the sepeiration point , the rate of vorticity shedding 
at the top is given hy the following equation, which replaces equation 




ds\ ^’u 
dt/ 5n 





(26) 


From equations (s) and (2l), 



Equations (27) yield the following expression for the right side of equa- 
tion (26): 



Expression (28) represents clockwise vortlcity shed at the top separa- 
tion point. At the hottom separation point, the amount of coimterclock- 
^se vorticity shed is also represented hy equation (28), if subscripts 
s are replaced hy s. 


Therefore, equating the net discharge of vorticity to zero. 


(uj)^ - (u?) 



2 J - 





f-i i)_ - («^ 3 


X29) 




In the quasi-steady flow, (u.^)— = and the coefficient of Ax must 

he zero. .Therefore, the coefficient of e in equation (29) must vanish: 


(sr - (sr 3^1 * ” f ■ f 35 ) 


= 0 
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or 


/ dq/^T] 





± - 

l^js 

[6r^da,J^ 


(R 

(R d’rr/6a)_ " ° 
s 


(30) 


Eq’iations (23), (25), and (30) and table I suffice to determine r (which 
appears in eq. (25)). The result is, 

r = -6.1 (31) 

Definition of separation point . - In effect, it has been assumed 
that, during the unsteady motion, separation is defined by the condition 
of zero shear (x = -0.1567) and the subsequent appearance of reverse 
flow relative to the surfane, Just as in steady flow. This assun^tion 
is open to question. The question is how (or whether) local velocity- 
profile characteristics may be interpreted to identify the leading edge 
of a wake. 


The usual steady criterion, which notes the appearance of reverse 
flow downstream of the point of zero shear, implies that the fluid in 
the wake is fixed relative to the body. If, in the unsteady case, the 
wake may still be regarded as fixed to the body, then it may be that the 
steady criterion is still applicable. 

However, the present assumption of the steady definition of separa- 
tion is not advanced with coniplete confidence. Rather, it is felt that 
only a suitable experiment can settle this point. 


Lift of Airfoil 

The steady lift coefficient of the airfoil of figure 1 is deter- 
mined from the steady circulation: 

q ( 0 ) ^ pur ^ r ^ 2 n( 0 . 076 l) = 0.48 ( 32 ) 

^ I PU2(2Z) ^ 


There are two contributions to lift proportional to e . From the un- 
steeidy circulation, 

=? 2rtr = -38.6 (33) 

*6 

and a further contribution is found by consideration of the remainder of 
the potential flow. The two components of Kelvin's impulse for the flow 
illustrated in figure 1 (leaving circulation out of account) are 

o ip 

1^,1 = jtpUZ (p cos a, sin a) 
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(See pars. 71 and 123 of ref. 11.) ‘Whence, the corresponding components 
of vector force are 

, Si 

= a — = jtpU2^a(-P^sin a, cos a) 


y 


and the lift is 


L = F^^cos a - F„sin a = jrpU2^a ^ ^ - 
y iL 


1 __ 

[ 1 + ^ cos 2a I 

1 + 3 "^ 


The following lift coefficient results: 


c(2)= I (1 + L ^ ^ ~ P 
^6 ^ I 1+3 


cos 2a I = 3.09 


(34) 


(34): 


The final expression for lift conibines equations (32), (33), and 


- =1°’ - 4h * "if) 

= 0.48 - 36 ^ 


(35) 


The sign of the second term of the result of equation (35) indicates 
that, while angle of atteick is increasing, the lift is lower than the 
qtiasi-steady value, and higher if the angle of attack is decreasing. 

Thus, near maximum lift, the lift curve would exhibit a counterclockwise 
hysteresis loop enclosing the stall point. This result is perhaps 
counter to expectations, because clockwise hysteresis is found experi- 
mentally for oscillating airfoils. It may be that different directions 
of hysteresis should be expected when the airfoil oscillates and when, 
aa in the present study, the stream direction oscillates. 


In any case, it may be shown that the overriding effect producing 
counterclockwise hysteresis in the present problem is the quasi-steady 
movement of the separation point over the top surface. As this separa- 
tion point moves forward under increasing angle of attack, clockwise 
vortlclty in the boundary layer joins the wake as the separation point 
passes. Accordingly, a counterclockwise airfoil circulation (negative 
y) is required to balance this effect. The term of equation (26) that 
is concerned with this movement is the one involving ds’/dt. 

For the elliptic airfoil problem treated herein, the quasi-steady 
movement of the upper separation point is quite extensive; numerically, 
^Tf/Sa = 13.9, indicating that the separation point is very loosely fixed 
to the airfoil surface. By way of contrast, in the case of the rotating 
circular cylinder, there is no effect of this sort, because a change in 
angle of attack produces no quasi-steady movement of the separation point 
at ail. Other contributions to the shedding of vortlcity then lead to 
the result of clockwise hysteresis for the circular cylinder. 
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LIFT HTSTERESIS 

The foregoing analysis does not provide a complete theory of un- 
steady flow about a stalled airfoil. Rather, the analysis illustrates 
the considerations that would underlie such a theory and further make 
plausible the general assurption of an expansion of lift coefficient in 
the form of equation (35). This expansion would be valid for nearly 
quasi-steady conditions. Also, it has been shown that a counterclockwise 
hysteresis locp may be expected at a nominal condition of maximum lift 
with, apparently, a large anplitude. 

If the foregoing conditions are met, and Ax is simple harmonic. 

Ax oc sin tot (36a) 


then 


xoc - 2 tod sin (2rtCJ0t - 90°) (36b) 

Use of relation (36b) in the lift formula (35) may be said to correspond 
to the assurption (ref. 5) of a positive lift-curve slope and a phase 
lag, 90° in this case, though the positive slope would not correspond to 
the steady lift curve. The hysteresis loop for this case appears as an 
ellipse on the lift curve of figure 6(a) . The amplitude and width of 
the loop are aasiamed small in the present discussion and are exaggerated 
in figure 6. 

If a is not sirple harmonic, then the concept of phase lag is al- 
together inappropriate. For example, if a is more nearly a "saw-tooth" 
function of time (fig. 6(b)), then the lift increment is nearly a 
"battlement" function. The corresponding hysteresis loop is nearly rec- 
tangular. If a changes according to an exponential pulse (illustrated 
by a Gaussian curve in fig. 6(c)), then the hysteresis loop is egg-shaped, 
with the broad end to the right. 


CQNCLUSIOKS 

The analysis of rotating stall in an axial-flow compressor requires 
specification of a dynamic lift curve applicable near stall. It has 
previously been suggested that unsteady lift hysteresis is an important 
charanteristlc of such a curve. 

Consideration of the familiar experimental fact of Magnus lift on 
a rotating cy3J.nder indicates a theoretical approach to the question of 
aerodynamiLc hysteresis which, though certainly not definitive, may prove 
helpful. The accepted explanation of Magnus lift is that, if the cylin- 
der is in motion toward the left and rotates clockwise, the movement of 


3778 



CF-S "back 8ll£ 


MCA TN 3571 


19 


the surfeice delays houadary-layer separation on the top and advances it 
on the hottom. In steady flow, considerations of constancy of wake cir- 
culation require that the outer velocity at the two separation points he 
equal. The delayed separation at the top inplies a lower velocity (vice 
versa on the hottom) , and a conpensatory clockwise circulation must 
therefore occiir. 

The foregoing reasoning is extended to apply to the problem of an 
airfoil of elliptic section in a stream of constant velocity hut of 
(slightly) oscillating direction. The airfoil is considered to he nom- 
inally at maximum lift. This assumption, reasonable for unsteady prob- 
lems at* nearly stalled conditions, provides an essential simplification. 
To first order in small quantities, the lift (circulation) increment due 
to the oscillation can depend only on rate of change of angle of attack; 
and, just as in the cylinder case, all induced wake effects may he 
ignored. 

For purposes of computing the amount of vorticity shed into the 
wake, the separation point is identified as the point of vanishing shear, 
just as in steady flow. It is not clear, however, that this assumption 
is proper. 

Under these various assumptions, the unsteady increment in lift 
coefficient of the ellipse is found to he -36 aZ/U, of which about 92 
percent is due to the unsteady movement of the separation points and the 
remiainder is due to impulsive pressure. This lift increment gives the 
amplitude of a lift -hysteresis loop at maxinnxm lift. The loop is coun- 
terclockwise, a result that can he related to the extremely migratory 
tendency of the separation point on the upper siorfeice of the ellipse 
under a change in angle of attaek in steady flow. 

Finally, assuming oscillations of low reduced frequency, certain 
observations may he made concerning the shapes of hysteresis loops, and 
the validity of the idea of a viscous time lag in connection with un- 
steady lift: If the angle of atteick undergoes harmonic oscillation, 

then the lift Increment is also harmonic with a 90° phase lead or lag, 

depending on the sign used in the definition of C 7 • In this case, the 

‘'e 

hysteresis loop is elliptic. If the angle of attack varies in a non- 
harmonic manner, then the variation of lift does not have the same de- 
pendence on time, and the idea of time lag is Inappropriate. 


Leid-s Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, -Ohio, August 17, 1955 



20 


KACA TN 3571 


APPENDIX - SYMBOLS 

Ttie following symbols are used in tMs report: 

A, a functions of ti (eqs. (?)) 

fiinctions of tj (eqs. (s)) 

C,c functions of tj (eqs. (l9b)) 

C, lift coefficient 

C, rate of change of lift coefficient with angle of attack a 

C-, . rate of change of lift coefficient with a 

‘'a, 

d function of t] (eq. (20b)) 

F universal function for houndary-layer calculation (eq. (5a)) 

f^ universal function for hovindary-layer calculation (eq. (l3)) 

2 semichord of elliptic cylinder (fig. l) 

n coordinate measured normal to surface (eqs. (s)^ (36)) 

q dimensionless outer velocity, 

R function of ti (eq. (3)) 

s dimensionless coordinate measured along surface of ellipse (fig. l) 

t time 

U stream velocity (a constant) 

u velocity parallel to surface 

W function related to unsteady houndary layer (eq. (l5)), = qZg 

X function related to unsteady houndary layer (eq. (is)) 

X Cartesian coordinate of surface of ellipse 

Y function related to unsteady houndary layer (eq. (l8)) 

y Cartesian coordinate of surface of ellipse 
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Z function related to momentum thickness of boundary layer (eq. (5a)), 

" VI 

a angle of attank 

a time rate of change of a 

A% Increment in angle of attack, = a, - 
P thickness ratio of elliptic cylinder 

r circulation in outer floa'j' 

TT coefficient of hysteresis in circulation (eq. (lO)) 

5 over-all thickness of boundary layer 

6* displacement thickness .of boundary layer 
e dimensionless angular velocity, = al/U 

T] coordinate on surface of ellipse (fig. 1, eq. (2)) 

8 momentum thickness of boundary layer 

X function for houndary-layer calculation (eq. (5a)) 

V kinematic viscosity coefficient 

p density 

T skin-friction coefficient 

<p velocity potential 

CD frequency 

Subscripts : 

s',£ evaluation at top or bottom separation point, respectively, of 
steady flow at maximum lift (a = ^ e = O) 

0 evaluation at forward stagnation point, where q = 0 

1 evaluation at outer edge of boundary layer 
evaluation at lower separation point 
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Subscript notation is used for partial differentiation where convenient 
Superscripts : 

(0) steady conditions at maximmn lif t 

(1) unsteady contribution due to movement of separation points 

(2) unstea^ contribution due to in^ulsive pressure 
evaluation at i:ipper separation point 

' denotes ordinary derivatives 
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TABLE I. - EESULTS AT SEPAEATIOW POINTS 


Quantity 

Separation point (x = -0.1567) 


Top (l?) 

Bottom (s) 

T1 

80.0° 

340.83° 

R 

0.985 

0.364 

q 

1.210 

-1.204 

S 

0.139 

0.137 

^Tja 

1.192 

-7.95 


-0.0062 

1.49 

Z 

1.11 

0.41 

a 

13.3 

-5.6 


-1.405 

0.52 

X 

-55.6 

-5.4 

Y 

4.9 

-4 . 2 



1'J 
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Figure 5. - Dietrlbution of peirameter Y of unsteady iDOundary layer. 
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Time, t 


Angle of attack, a 


Figure 6. 


(a) Harmonic a(t). 

- Variation of lift coefficient and angle of attack with time and with each other. 
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Tirnsj t 


An^e of attack, a 


Figure 6. 


(b) "Saw-tooth" function a(t). 

- Contlmied. Variation of lift coefficient and angle of attack with time and with each other. 
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(c) Exponential pulse 


Figure 6 


- Concluded 


Yeoriatlon of lift coefficient and 


coefficient 



of attack with time and with each other 



